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Abstract. Selected topics in solar dynamo theory are being highlighted. The possible relevance 
of the near-surface shear layer is discussed. The role of turbulent downward pumping is men- 
tioned in connection with earlier concerns that a dynamo-generated magnetic field would be 
rapidly lost from the convection zone by magnetic buoyancy. It is argued that shear-mediated 
small-scale magnetic helicity fluxes are responsible for the success of some of the recent large- 
scale dynamo simulations. These fluxes help in disposing of excess small-scale magnetic helicity. 
This small-scale magnetic helicity, in turn, is generated in response to the production of an 
overall tilt in each Parker loop. Some preliminary calculations of this helicity flux are presented 
for a system with uniform shear. In the Sun the effects of magnetic helicity fluxes may be seen 
in coronal mass ejections shedding large amounts of magnetic helicity. 
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■ 1. Introduction 

Unlike the geodynamo, which is widely regarded a solved problem (e.g. Glatzmaier & 
Roberts 1995), the solar dynamo problem remains unsolved in that there is no single 
model that is free of theoretical shortcomings and that actually reproduces the Sun. The 
flux-transport dynamo models (e.g. Dikpati & Charbonneau 1999; Kiiker et al. 2001; 
Chatterjee et al. 2004) are currently the only models that display some degree of realism, 
but there are aspects in their design that are arguably not sufficiently plausible. One 
problem is the assumption of a rather low turbulent magnetic diffusivity that needs to 
be assumed in a ad hoc fashion. On theoretical grounds the turbulent magnetic diffusivity 
should be comparable with the turbulent kinematic viscosity (Yousef et al. 2003), which 
Q\ | is not the case on these flux-transport dynamo models. 

The purpose of this review is to collect recent findings that are relevant in solving 
the solar dynamo problem. On the one hand, there are direct simulations of convective 
dynamo action in spherical shells by Brun & Toomre (2002), Brun et al. (2004, 2006), 
Browning et al. (2006), and Brown et al. (2007). These models are already quite realistic 
and should ultimately be able to reproduce the solar dynamo. Their main shortcom- 
ings range from still insufficient resolution (even though it is already exhausting current 
capabilities) to the negligence of certain physical features of the model. Examples of 
shortcomings include the still insufficient degree of stratification as well as the simpli- 
fied treatment or even the neglect of surface and tachocline boundary layers. On the 
other extreme, there are phenomenologically motivated mean-field models that are con- 
structed based on their success in reproducing the Sun. These models tend to utilize a 
subset of known turbulent transport effects with ad hoc amplitudes and prescriptions of 
their radial profiles. Such models arc made nonlinear by simple a quenching terms. How- 
ever, this type of prescription that is still commonly used in mean field models does not 
correctly describe the saturation behavior of large-scale magnetic fields as known from 
three-dimensional turbulence simulations. This is true even for much simpler models, 
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Figure 1. Comparison of the differential rotation contours that were originally expected by 
Yoshimura (1975) based on solar dynamo model considerations (left) with those by Thompson 
et al. (2003) using helioseismology (middle). Note the similarities between the contours on the 
left (over the bulk of the convection zone) and those on the right (over the outer 5% of the 
solar radius). In the right hand panel we show radial profiles of angular velocity as given by 
Benevolenskaya et al. (1999). Note the sharp negative radial gradient near the surface. 

for example those where the turbulence is driven at a specific length scale via a forcing 
function (Brandenburg 2001). 

2. Negative radial shear? 

In the 1970s a number of solar dynamo models were discussed that were based on an 
a effect that is positive in the northern hemisphere (and negative in the southern hemi- 
sphere) and a differential rotation that had a negative radial gradient, i.e. dil/dr < 
(see the left panel of Fig. [T|) . Such models produce an equatorward migration of mag- 
netic activity toward the equator. Examples of such models in full spherical geometry 
include the papers Steenbeck & Krause (1969), Roberts & Stix (1972), Kohler (1973), 
and Yoshimura (1975). During the 1980s these models faced two difficulties, a conceptual 
one and one based on a conflict with observations. The conceptual difficulty has to do 
with the idea that the magnetic field may be fibral (Parker 1984), i.e. it is distributed in 
the form of many flux tubes with a small filling factor and relatively large field strength. 
However, such fields would easily be unstable to magnetic buoyancy instabilities. This 
and other considerations led Spiegel & Weiss (1980) to propose that the solar dynamo 
might instead work at the bottom of the solar convection zone and not, as assumed until 
them, in a distributed fashion in the bulk of the convection zone. In fact, in an earlier 
paper by Parker (1975) this possibility was already discussed in some detail. The conflict 
with observations had to do with the fact that helioseismology began to put tight con- 
straints on the form of the Sun's internal angular velocity profile. This seemed to rule 
out a negative radial Q gradient. 

These two difficulties led to a new type of dynamo model that works in the lower 
overshoot layer, where a has the opposite sign (DeLuca & Oilman 1986, 1988). Not 
only would magnetic flux tubes presumably be stabilized against magnetic buoyancy 
instabilities, but at that depth the sign of the a effect would be reversed and thus dynamo 
waves would again propagate equatorward even though dfl/dr > 0. However, several 
problems with this approach have been discussed (phase relation, number of toroidal 
flux belts, etc.). Yet another model is the flux transport model, where a meridional 
circulation (equatorward at the bottom of the convection zone) is chiefly responsible for 
driving the dynamo wave equatorward. However, it is now quite clear that there is a 
strong shear layer near the surface where dfl/dr < 0, just as originally anticipated. The 
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Figure 2. Horizontally averaged magnetic fields B x (upper panel on the left) and B y (lower 
panel on the left), as functions of time and z for a run with shear and no rotation (Kapyla et 
al. 2008), as well as B y from a local model of an accretion disc (right hand panel), where z = 
corresponds to the midplane. On the left-hand panel the dotted white lines show top (z — d) 
and bottom (z = 0) of the convection zone. 

difference is that now this shear layer only extends over the outer 5% of the solar radius. 
In the past the relevance of this near-surface shear layer was discarded mainly on the 
grounds that magnetic buoyancy effects would lead to a rapid loss of magnetic field. 
Another reason is that near the surface the local turnover time is still short compared 
with the rotation period. However, simulations of compressible convection did show some 
time ago that magnetic buoyancy is effectively overpowered by the effects of turbulent 
downward pumping. Regarding the timescales, rotational effects are likely to play a role 
at the bottom of the near-surface shear layer, which is at a depth of about 40 Mm. The 
relative importance of rotation and convection is determined by the Coriolis number, 

Co = 2flr, (2.1) 

where fl is the solar rotation rate and r = H p /u lms is the turnover time with H p being the 
pressure scale height. Using il = 3 x 10 _6 s _1 for the solar rotation rate, w rms = 50m/s 
for the rms velocity at a depth of 40 Mm, and H p = 13 Mm for the pressure scale height 
at that depth we find r = 1.3 d, and hence Co = 1.4, suggesting that rotation should 
become important at the bottom of the near-surface shear layer. 




3. Quenching and examples of large-scale dynamos 

Since the early 1990s there has been the notion that catastrophic a quenching has 
been a major concern in dynamo theory (Vainshtein & Cattaneo 1992). This led Parker 
(1993) to suggest that the solar dynamo might work in a segregated fashion where shear 
acts only near the lower overshoot layer where the turbulent magnetic diffusivity is low. 
In the bulk of the convection zone, on the other hand, the magnetic field is supposed to 
be weak, and so there was hope that the a effect would not be catastrophically quenched. 
However, only later the origin of catastrophic quenching was understood to be due to 
magnetic helicity conservation. As a consequence, catastrophic quenching could not be 
alleviated by spatial rearrangements, but only by magnetic helicity fluxes out of the 
domain. Indeed, there are now a number of simulations that show successful generation 
of large-scale magnetic fields; see Fig. [2] All these simulations have in common that 
there are contours of constant shear crossing the outer surface. This is believed to be 
important because the magnetic helicity flux is expected to be directed along the contours 
of constant shear (Brandenburg & Subramanian 2005a). 
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Figure 3. Original sketch of a Parker loop (Parker 1957) on the left, compared with a loop 
constructed from the Cauchy solution of Yousef & Brandenburg (2003) and its corresponding 
magnetic helicity spectrum Ht on the right. 

4. a effect and production of small-scale magnetic helicity 

The a effect is central to many mean-field approaches of the solar dynamo. Its pres- 
ence is usually associated with the systematic action of the Coriolis force on convecting 
fluid elements. As such an element rises vertically, it expands and counter rotates due to 
the Coriolis force, thereby attaining a downward oriented vorticity in the northern hemi- 
sphere, and an upward oriented vorticity in the southern hemisphere. As the magnetic 
field is dragged along with the gas, a systematic poloidal field component is being pro- 
duced from a toroidal field. An early sketch of this process was given by Parker (1957), 
where a rising loop attains a systematic tilt; see the left hand panel of Fig. [3] 

In recent years the nonlinear saturation of this process has been understood to be due 
to the internal twist that is being induced as the loop begins to tilt. This is already quite 
evident from the early sketch by Parker (1957) suggesting the presence of a spatially 
extended structure of the loop with oppositely oriented twisting motions on both ends 
of the loop; see the round arrows on the tube in the left hand panel of Fig. [3l 

In the middle panel of Fig. [3] we show a visualization of an analytically generated rep- 
resentation of a Parker loop using the Cauchy solution for ideal magnetohydrodynam- 
ics together with the corresponding numerically constructed magnetic helicity spectrum 
(Yousef & Brandenburg 2003); see the right hand panel of Fig. [3j This spectrum shows 
quite clearly a negative magnetic helicity at low wavenumbers k/k\ ~ 1 (where k\ = 2n / L 
is the smallest wavenumber in a box of scale L) and positive magnetic helicity at larger 
wavenumbers around k/k\ ~ 4. This confirms an early finding of Seehafer (1996) that 
the a effect corresponds to a segregation of magnetic helicity in scale and that magnetic 
helicity of smaller scale and opposite sign is inevitably being produced by the a effect. It 
is this what quenches the a effect in a potentially catastrophic manner. Alleviating this 
type of quenching requires that we have to get rid of this smaller-scale magnetic helicity 
through helicity fluxes. This is something that the Sun seems to accomplish through 
coronal mass ejections (Demoulin et al. 2002, Blackman & Brandenburg 2003). In the 
following we describe in more detail the magnetic helicity flux and its relation to the 
amount of shear. 

5. Magnetic helicity flux 

There are a number of analytic calculations of magnetic helicity fluxes relevant to 
astrophysical dynamos (Kleeorin et al. 2000, 2002, 2003a,b; Vishniac & Cho 2001; Sub- 
ramanian & Brandenburg 2004, 2006; Brandenburg & Subramanian 2005a, b). Here we 
present preliminary numerical calculations in a simple system. We calculate the small- 
scale current helicity flux, IFq r, in homogeneous turbulence in the presence of shear, S, 
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Figure 4. B dependence of — CVc f° r -Bo = Bo y y with i? m = 1.4, Sh pa —2 (left), 5 dependence 
of Cvc for _R m = 30, Bo/B cq ~ 0.2 (middle), and i? m dependence of Cvc for Sh ~ 0.5, 
Bo/B cq pa 0.2. 

and a uniform magnetic field _B - Since the system is completely uniform, it makes sense 
to consider full volume averages, denoted here by an overbar. The wavenumber of the 
averaged quantities is zero, so we have infinite scale separation and can then force the 
turbulence at the scale of the system, i.e. at wavenumbers between 1 and 2 times the 
smallest finite wavenumber, k\. Our average forcing wavenumber is therefore fcf/fci = 1.5. 
We consider the dependence on the parameters Bq/B qci , i? m and Sh, where 

B cq = Mo(p« 2 }, Rm = u rms /^kf, Sh = S/i]kf. (5.1) 

The small-scale current helicity flux is calculated from the expression (Brandenburg & 
Subramanian 2005b) 

IFl? = 2eVj + (V x e) x b, (5.2) 

where 

e = -u x (b + B ) + r^ j (5.3) 

is the small-scale electric field and j = V x b/fj,o is the small-scale current density with 
Ho being the vacuum permeability. In order to avoid taking more than two derivatives 
we integrate the second expression in equation (|5.2|) by parts, i.e. 

(V x e) x b = (V6) T e, (5.4) 

where we have used the fact that V • b — 0. The i component of this term can also be 
written as bj^ej, where a comma denotes partial differentiation. 

Based on calculations using the minimal r approximation the small-scale current he- 
licity flux is expected to be given by (Brandenburg & Subramanian 2005a) 

= Cvc (SB) xB, (5.5) 

where Cvc is a non-dimensional coefficient that is of the order of St 2 , where St = Tu Trns k{ 
is the Strouhal number. Throughout this work we use a uniform shear flow, i.e. U = 
(0, Sx, 0), so the z component of J-q is 

y S c z = \C vc s(B 2 y -B 2 x ). (5.6) 

In the following we quote values of Cvc that are computed by imposing a uniform field 
either in the x or in the y direction. The numerical resolution is only 32 3 . In the left hand 
panel of Fig. |4] we show the dependence on B for a small value of i? m (i? m = 1.4) and 
fixed shear parameter Sh pa —2 (for S < 0). So far there is no indication that the flux 
depends on Sh; see the middle panel of Fig. [4] However, the flux shows a strong increase 
with R m ; see the right hand panel of Fig. |4j It should be noted that, since we keep the 
forcing unchanged while changing the viscosity, the resulting rms velocity also changes, 
and so the resulting values of Sh and Bo/B cq also change somewhat. 
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6. Discussion 

In this paper we have presented selected aspects of the solar dynamo problem where 
there has been recent progress. The idea that the solar dynamo may operate in the bulk 
of the convection zone is motivated in part by the fact that it is easier to dispose of 
small-scale magnetic helicity from upper layers than from deeper down. Also the sign 
of the radial differential rotation is negative and would produce equatorward migration 
of dynamo waves in the presence of a positive a effect in the northern hemisphere of 
the Sun. A number of other arguments for a distributed solar dynamo, where the field 
that makes sunspots does not solely come from the lower overshoot layer, have been 
discussed elsewhere (Brandenburg 2005). However, what is not yet well addressed is 
the production of sunspots within active regions, and perhaps even the active regions 
themselves. It is plausible, and it has indeed been argued, that sunspots can be the 
result of negative turbulent magnetic pressure effects (Kleeorin & Rogachcvskii 1994) or 
a turbulent magnetic collapse phenomenon (Kitchatinov & Mazur 2000). Both processes 
rely on turbulent transport processes that could be verified numerically and whose effects 
could also be demonstrated directly in suitably arranged simulations. 
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Discussion 

Strassmeier: Would a small (ideal) isolated spot trace the large-scale or the small-scale 
field? 

Brandenburg: Large-scale field means really mean field, and if we define this through 
azimuthal averaging, not much of the spot field would survive azimuthal averaging. There- 
fore, most of the spot field would contribute to the small-scale field and only a small 
fraction would contribute to the large-scale field. 



